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Summary. We obtain necessary optimality conditions for a semi-discretized opti- 
mal control problem for the classical system of nonlinear partial differential equations 
modelling the water-oil (isothermal dead-oil model). 
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1 Introduction 

We study an optimal control problem in the discrete case whose control system 
is given by the following system of nonlinear partial differential equations, 



which result from a well established model for oil engineering within the frame- 
work of the mechanics of a continuous medium [3] . The domain I? is an open 
bounded set in with a sufficiently smooth boundary. Further hypotheses 
on the data of the problem will be specified later. 

At the time of the first run of a layer, the flow of the crude oil towards 
the surface is due to the energy stored in the gases under pressure in the 
natural hydraulic system. To mitigate the consecutive decline of production 
and the decomposition of the site, water injections are carried out, well before 
the normal exhaustion of the layer. The water is injected through wells with 
high pressure, by pumps specially drilled to this end. The pumps allow the dis- 
placement of the crude oil towards the wells of production. More precisely, the 
problem consists in seeking the admissible control parameters which minimize 



dtu — Alp(u) = div (g(u)Vp) 



in = i7 X (0,T), 
inQr^nx (0,T), 



dtp - div ((i(M)Vp) = / 



(1) 
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a certain objective functional. In our problem, the main goal is to distribute 
properly the wells in order to have the best extraction of the hydrocarbons. 
For this reason, we consider a cost functional containing different parameters 
arising in the process. To address the optimal control problem, we use the 
Lagrangian method to derive an optimality system: from the cost function we 
introduce a Lagrangian; then, we calculate the Gateaux derivative of the La- 
grangian with respect to its variables. This technique was used, in particular, 
by A. Masserey et al. for electromagnetic models of induction heating [1, 7], 
and by H.-C. Lee and T. Shilkin for the thermistor problem [5]. 
We consider the following cost functional: 

j{u,p, f) = \ \\u - u\\Iq^ + \\\p- p\g^Q^ + 1 wfwii^Q^ + f mwiQ^ . 

(2) 

The control parameters are the reduced saturation of oil u, the pressure p, 
and /. The coefficients /3i > and /32 > are two coefficients of penaliza- 
tion, and go > 1- The first two terms in (2) allow to minimize the difference 
between the reduced saturation of oil u, the global pressure p and the given 
data U and P. The third and fourth terms are used to improve the quality of 
exploitation of the crude oil. We take ^2 = just for the sake of simplicity. It 
is important to emphasize that our choice of the cost function is not unique. 
One can always add additional terms of penalization to take into account 
other properties which one may wish to control. Recently, we proved in [8] 
results of existence, uniqueness, and regularity of the optimal solutions to the 
problem of minimizing (2) subject to (1), using the theory of parabolic prob- 
lems [4, 6]. Here, our goal is to obtain necessary optimality conditions which 
may be easily implemented on a computer. More precisely, we address the 
problem of obtaining necessary optimality conditions for the semi-discretized 
time problem. 

In order to be able to solve problem (l)-(2) numerically, we use discretiza- 
tion of the problem in time by a method of finite differences. For a fixed real 
A'', let T = be the step of a uniform partition of the interval [0, T] and 
tn = riT, n= 1,. . . ,N. We denote by an approximation of u. The discrete 
cost functional is then defined as follows: 

N 

j{u\p\n = \Y. / {\w''-nln + \\p^-nln+Mr\\Z,n)dx. 

(3) 

It is now possible to state our optimal control problem: find /") which 

minimizes (3) among all functions {u",p", f") satisfying 

' - = div {g{u'')Vp) 

^ P''^'-P" - div (d(u")Vp") = /" 

.Plaf2 = 0' P\t=o= Po- 



rn fl , 
in i? , 



(4) 
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The soughtafter necessary optimality conditions are proved in §3 under suit- 
able hypotheses on the data of the problem. 

2 Notation, hypotheses, and functional spaces 

Our main objective is to obtain necessary conditions for a triple /") 
to minimize (3) among all the functions /") verifying (4). In the sequel 

we assume that ^p, g and d are real valued functions, respectively of class C^, 
and C^, satisfying: 

(HI) < ci < d{r), ^(r) < ca; \d'{r)\, \ip'{r)\, |<^"(r)| < C3 Vr G K. 

(H2) uq, po G (n), and U, P e L'^{0), where uq, po,U, P : Q ^ M, and 

MdQ = Polaa = 0- 
We consider the following spaces: 

w^{n) := {u G u>in), vu G L^in)} , 

endowed with the norm \\u\\^^f^^^^ = \\u\\j,jj + ||Vu||p_^^; 

W^{f2) := {u G w^{n), V\ G LP{n)} , 
with the norm ||w||v^'2(/3) = ||^i|lvKi(fi) + ll^^^llp r?' ^^'^ following notation: 

W :=W2g {O) ; 

T := L29(|2) ; 

2-i 

H:=L'">{n)x W2g" {^). 

3 Main results 

We define the following nonlinear operator corresponding to (4): 

F -.W xW xY — > H xH 
iu\p\n^Fiu\p\n, 

where 



F{u^,p^,r) = 




^ - Z\(/?(m") - dw(g(u")Vp"), 7ou" -uo\ 
^ - div {d{u-)Vp'-) - r, -Po) ' 



7o being the trace operator 7om" = w|t=o- hypotheses ensure that F is 
well defined. 
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3.1 Gateaux differentiability 

Theorem 1. In addition to the hypotheses (HI) and (H2), let us suppose that 
(H3) |^"'| < c. 

Then, the operator F is Gateaux differentiable and for all (e, w,h) G W x 
W xY its derivative is given by 

iF(o",p", /»)(e, w,h) = ^F (»" + ae, p" + »«., /" + »/.) 

= (*F.*F.)^(|;|). 

^1 = e-dw(¥''(w")Ve)-dw(¥'"(w")eVw") -div (3(w")Vu;) -div (5'(w")eVp"), 
^2 = 7oe, ^3 = w — div {d{u")Vw) — div (d'(u")eVp") — h, ^4, = 70W. 
Furthermore, for any optimal solution (-(1",^",/") of the problem of mini- 
mizing (3) among all the functions {u",p"',f"') satisfying (4), the image of 
6F /") is equal toHxH. 

To prove Theorem 1 we make use of the foUowing lemma. 

Lemma 1. The operator 6F{u'^,p'^, f"^) : W x W x T — > H x H is linear 
and bounded. 

Proof (Lemma 1). For all {e,w,h) gW xW xT 

^„nFi(u",p",r)(e,«;,/i) 

= e-div ((^'(M")Ve) - div (y^"(M")eVu") 
- div {g{u''-)Vw) - div (ff'(tt")eVj9") 
= e - (p'(u")Ae - (^"(u")Vu".Ve - (^"(u")eAM" 
- (^"(w")Ve.Vu" - <^"'(u")e|Vu"|^ - £f(u")Au; - 5'(u")Vu".Vw; 

- 5'(u")eAp" - g'(u")Ve.Vp" - 5"(u")eVu".Vp" , 

where ^^ni^" is the Gateaux derivative of F with respect to m". Using our 
hypotheses we have 

\\g"{vP)eVvr.Vp^\\2,,n < \\e\\oo,n\\^vP ■^P'^Wn 

<\\eU,n\\^u-\\^,a\m\A,a 

< c||u"||w'||p"|kl|e||iv. 
Evaluating each term of 5„nFi, we obtain 



\\Su^F,{u\p\f"){e,w,h)hg,Q^ 

<c{\\u-\\w, IKIk, lirilr) {\\e\\w + Mw + Mr) ■ (5) 
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In a similar way, we have for all {e,w,h) G W x W x T that 

(5pnF2(w",p",/")(e,w,/i) = w-div (d{u'')Vw) - div {d' {u")eVp'') - h 
= w- d{u")Aw - d'{u")Vu'\\7w - d'(M")eAp" 

- d'(M")Ve.VM" - d'(?i")eVu".Vp" - h , 

with SpnF the Gateaux derivative of F with respect to p". Then, using again 
our hypotheses, we obtain that 

\\Sp^F2{u'\p", Die, w, h)\\2g,n < \\w\\2g,n + \\Vw\\2g,a + c\\Aw\\2g,a 
+ c\\Vu^.Vw\\2g,o + c\\eAp^2g,n 
+ c\\Ve.Wu^2g,n + c||eVu".Vp"||2g,fi + \\h\\2g,n ■ (6) 

Applying similar arguments to all terms of (6), we then have 

\\5p.F2{u^,p^,n{e,w,h)\\2,,f2 

<c{\\u-\\w, \\p-\\w, lirilr) (||e|k + Mw + Mr) ■ (7) 

Consequently, by (5) and (7) we can write 

\\5F{u^,p^,n{e,w,h)\\H>cH.r 

< c(||u"|k, IKIk, \\n\r) {\\e\\w + Mw + \\h\\r) • 

□ 

Proof (Theorem 1). In order to show that the image of dF{u,p, f) is equal to 
H X H, we need to prove that there exists {ejWjh) € W x W x T such that 

e - div {ip'{vP)Ve) - div {if" {W)e\/TF) 

— div (g(u")Vw) — div (g'(u")eVp") = a , 
w - div (d(u«)Vu;) - div (d'(M")eVp") -h = P, (8) 



e 



on 



0, e\t^Q = b, 



for any (a, a) and (J3, b) € H . Writing the system (8) for /i = as 

e - (p'(u^)Ae - 2(^"(it^)Vlt^.Ve - ip'' (^)eAW - ip''' (^)e\S/vJ^\^ , 
-g{W)Aw - g'{v^)VW.Vw - g'{TF)eAp^ 

- c/'(u")Vp"-Ve - 5"(tI")eVu".Vp" = a , 
w - d{vP)Aw - d'(vP^)'SlW.Vw - d'(vP)eAp^ (9) 
- VtI".Ve - d'(tI")eVu".Vp" = (3 , 
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it follows from the regularity of the optimal solution that <f"{u")Au"', 
(^"'(u")|Vit"p, g'{T^)Ap^, .g"(it")Vit".Vp", d'(u")Ap", and d' (v^)Vv^ .Vjf^ 
belong to L2«o(I2); (^"(it»)Vit", g'(i!^)Vv^, g'{v/^)\/p^, and d'(M")VM» be- 
long to L'^'^°{n). This ensures, in view of the results of [4, 6], existence of a 
unique solution of the system (9). Hence, there exists a {e,w,0) verifying (8). 
We conclude that the image of SF is equal to H x H. □ 

3.2 Necessary optimality condition 

We consider the cost functional J : W x W x T ^ R {3} and the Lagrangian 
C defined by 

C iu^,p^, r,puei, a,b) = J {u^,p-, /") + (^F{u",p-, /"), 

where the bracket (•, ■) denotes the duality between H and H'. 

Theorem 2. Under hypotheses (H1)-(H3), if(}^,W^f") is an optimal solu- 
tion to the problem of minimizing (3) subject to (4), then there exist functions 
(ei,pi) € (J7) X W2{f^) satisfying the following conditions: 

el + div ((/9'(u") Vei) - d'(w")Vp".Vpr - iy9"(u") Vu". VeT 

JV 

- fi('(tI")Vp".Ver = r ^(u" - U) , 

n=l 

ei\en = , e^l^T = ' 

_ _ ^_ (10) 

+ div ((i(M")Vpr) + div {g{u^)^el) = t - P) , 

n=l 

Pilafi^O' Prlt=T = 0' 

JV 

9o/3iT^ir|2«°-27^ = pr. 

Proof. Let [u",p", f") be an optimal solution to the problem of minimizing 
(3) subject to (4). It is well known (cf. e.g. [2]) that there exist Lagrange 
multipliers ({pi,a), {eT,b)) G H' x H' verifying 

^(„„,p„j„)/:(i^,^,7^,pr,er,o,6) ie,w,h) = y{e,w,h) eWxWxY, 

with (5(u'',p",/")/^ the Gateaux derivative of C with respect to (u",p",/"). 
This leads to the following system: 
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N f _ _ 

rV / ((^-f/)e + (^-P)u; + go/3i|/"P*"V"/i) (^a; 

- j i^e - div (v?'(M")Ve) - div (<^"(u")eVM") 

— div {g(v^)Vw) — div ((?'('u")eVp")^ eT^ dx 

— I {w — div (d(u")Vti;) — div {d' (vP)eSIp^) —h)pidx 

J Q 

-(7oe,a) + -(7ow,5) = \/{e,w,h) £W xW xT. 
The above system is equivalent to the following one: 

/ I T J7)e - ((i'(u^)eVp^)pr+ eel - div (<^'(it^)Ve) el 
-div (v"(ti")eVM") eT - rfw (5'(ti")eVp") ej^ da; 
+ / I r /^(p" — -P)w + wpT — div (d(u")Vw) pT — dw {g(vP)Vw) eT ] da; 

•^^ V n=l / 

+ (^go/3iTXj|7^|2«°-27^/i-pr/ij da; 
+ (7oe,a) + (70 w, 6) = V(e,w,/i) G x x T. 

(11) 

In others words, we have 

/ I r - J7) + d'(u)Vp^.Vpr - eT - dw ((y9'(it^)VeT) 

+ (p"(u^)Vu^.Ver + g\u^)Vp^SJelj e da; 

r ( ^ \ (12) 

+ y f r ^(p^ - P) + pT - diw (d(z?") VpT) - div (g (m^) Vel) j w; da; 

+ |^go/3iT X] iW^-^rh - pThj dx 

+ (7oe,a) + (70 w, 6) = y{e,w,h) eW xW xT 
Consider now the system 



Moulay Rchid Sidi Ammi and Delfim F. M. Torres 

ei + div {ip'{v/')Vei) - rf'(w")Vp"-Vpi - (/3"(tI")Vti".Vei 



N 



- 5('(u")Vp".Vei = T ^(w" - U) , 



n=l 

AT 



(13) 



Pi + div (d(M")Vpi) + div (g(M")Vei) = t ^(p" - P), 



n=l 



eiUfi = Pilafi = 0> eilj^y = pilj^j, = 0, 

with unknowns (ei,pi) which is uniquely solvable in W|(i?) x W|(/2) by 
the theory of elliptic equations [4]. The problem of finding (ejW) G W x W 
satisfying 

e — div ((/?'(M")Ve) — div ((^"(it")eVu") — div (,g(u")Vw) 

— div ((7'(u")eVp") — sign{ei — e^) , 

_ (14) 

w — div {d{u")'Vvu) — div {d (zt")eVp") = sign{pi — pi) , 

Joe = jqw = , 

is also uniquely solvable on W|^(I2) x W|^(i7). Let us choose /i = in (12) and 
multiply (13) by {e,w). Then, integrating by parts and making the difi^erence 
with (12) we obtain: 

y" (^e - div (<^'(ii")Ve) - div (</9"(u")eVw") - div {g{vF)'^w) 

— div ((7'(M")eVp")^ (ei — el) dx ^^^^ 

+ {w — div ((i(u")Vw) — div (d'(ti")eVp")) (pi — pi) dx 
Jo 

+(7oe, 7oer - a) + {jqw, ^oPi -b) = V(e, w) gW xW. 
Choosing (e, tw) in (15) as the solution of the system (14), we have 

/ sign{ei — er)(ei — eT) dxdt + / sign{pi — pi){pi — pi) dx = . 
Jn Jn 

It follows that ei = eT and pi = pi. Coming back to (15), we obtain 7061 = a 
and 7opr = b. On the other hand, choosing (e, w) = (0, 0) in (12), we get 

/ (^PiTj2j7^\'^''°~^7^-p?jhdx = 0,Wher. 

Then (10) follows, which concludes the proof of Theorem 2. □ 

We claim that the results we obtain here are useful for numerical im- 
plementations. This is still under investigation and will be addressed in a 
forthcoming publication. 
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